library(rstatix)
library(dplyr)
library(ggplot2)
# A small function helping us to round all numbers (0.123456 -> 0.123)
round_df <- function(df, digits = 3) {
nums <- vapply(df, is.numeric, FUN.VALUE = logical(1))
df[,nums] <- round(df[,nums], digits = digits)
(df)
}
# Loading the CSV file
df <- read.csv("data-wpm-threefactorial.csv")
# Show header if everything is OK
head(df)
# Check if data is normal distributed
round_df(df %>% group_by(Swipe, Visual, Tactile) %>% shapiro_test(WPM))
# Nothing is significant, so we can run a three-way repeated-measures ANOVA
anova <- anova_test(data = df, dv = WPM, wid = SubjectID, within = c(Swipe, Visual, Tactile), effect.size = "pes")
# Auto correction of degree of freedom (if necessary)
get_anova_table(anova, correction = "auto")
# Let's aggregate the data
means <- df %>% group_by(Swipe,Tactile, Visual) %>% summarize(
Mean = mean(WPM),
SD = sd(WPM),
N = n(), .groups = 'drop',) %>% mutate(
SE = SD / sqrt(N),
CI = qt(1 - (0.05 / 2), N - 1) * SE)
# Let's plot the data
plot <- ggplot(means, aes(y=Mean, x=Visual, fill=Swipe)) +
geom_bar(stat = "identity", position = position_dodge(width = 0.8), width=0.75) +
geom_errorbar(aes(ymin=Mean-CI, ymax=Mean+CI, width=.1), position=position_dodge(0.8), width=0.2, size=0.3) +
ylab("words per minute") +
xlab("")+
ggtitle("Typing performance in words per minute (WPM)") +
theme_minimal() +
facet_wrap(Tactile~.) +
scale_fill_manual(values = c("tomato1","steelblue2","gold2", "olivedrab3", "gray70")) +
theme(
legend.margin=margin(t = -5, unit='mm'),
plot.margin=unit(x=c(1,2,1,1),units="mm"),
legend.position = "bottom",
legend.title = element_blank(),
plot.title = element_text(hjust = 0.5),
panel.background = element_rect(size=0.5, color="transparent", fill="gray98" )
)
options(repr.plot.width = 6, repr.plot.height = 4)
plot
ggsave("Swipe-Modality-BarChart.pdf", plot, width=5, height=3.5, device=cairo_pdf)
